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Abstract 

We present a search for the decays t~V and K~vV in a 253 fb~^ data sample 

cohected at the T(4S') resonance with the Belle detector at the KEKB asymmetric energy B factory. 
Combinatorial and continuum backgrounds are suppressed by selecting a sample of events with one 
fully reconstructed B. The decay products of the B on the other side of the event are analyzed to 
search for B~ t~V and B~ — > K~vV decays. We find no significant evidence for a signal and set 
{) confidence level upper limits of B{B- t'V) < 1.8 x 10"^ and B{B- R-vV) < 3.6 x 10"^ 



All results are preliminary. 
PACS numbers: 
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The purely leptonic decay B~ (charge conjugate states are imphed throughout 

the paper) is of particular interest since it provides direct measurement of the product of the 
Cabibbo-Kobayashi-Maskawa(CKM) matrix element Vub and the B meson decay constant 
fs- In the Standard Model(SM), the branching fraction of the decay B~ i^V is given as 

Stt y J 

where is the Fermi coupling constant, mi and rriB are the charged lepton and B meson 
masses, tb is the B~ lifetime. The dependence on the lepton mass arises from helicity 
conservation, which suppresses the muon and electron channels. The CKMfitter predicts 
the B~ r~z/ branching fraction to be (9.3 tiJ) x 1O~^0. No evidence for an enhancement 
relative to the SM prediction was observed in previous experimental studies. The most 
stringent upper limit has been achieved by the BABAR Collaboration : B{B~ t~V) < 
4.2 X 10"^ at 90% confidence level (C.L.) 0. 

Flavor-changing neutral-current transition such as 6 — > sui? occurs in the SM via one- 
loop box or electroweak penguin diagrams with heavy particles in the loops. Because heavy 
non-SM particles could contribute additional loop diagrams, various new physics scenarios 
can lead to significant enhancements in the observed rates MJ^I- The SM B~ — > K'uu 
branching fraction has been estimated to be (S.Slgie) ^ 10~^|5|, |6[, while the most stringent 
pubhshed experimental limit is B{B~ R-uV) < 5.2 x 10"^ at 90% C.L. 

We use a 253 fb~^ data sample containing 275 x 10^ B meson pairs collected with the 
Belle detector at the KEKB asymmetric energy e~^e~ (3.5 on 8 GeV) collider operating 
at the T(4S') resonance {^/s = 10.58 GeV). The Belle detector is a large-solid-angle mag- 
netic spectrometer consisting of a three-layer silicon vertex detector, a 50-layer central drift 
chamber (CDC), a system of aerogel threshold Cerenkov counters (ACC), time-of- flight scin- 
tillation counters (TOF), and an electromagnetic calorimeter comprised of Csl(Tl) crystals 
(ECL) located inside a superconducting solenoid coil that provides a 1.5 T magnetic field. 
An iron fiux-return located outside of the coil is instrumented to identify K^^ and muons. 
The detector is described in detail elsewhere j^. 

The strategy adopted for this analysis is to reconstruct exclusively the decay of one of 
the B mesons in the event and compare properties of the remaining particle (s) in the event 
(referred to as the signal side) to those expected for signal and background. All the tracks 
and photon candidates in the event not used to reconstruct the B are studied to search for 
B^ — s> T^V and B~ K~vV. The advantage of having a sample of fully reconstructed B 
meson is to provide a strong suppression of the combinatorial and continuum background 
events. The disadvantage is the low efficiency of full B meson reconstruction (approximately 
0.3%). 

Fully reconstructed B mesons, -Brec? are observed in the following decay modes: B^ — >■ 

;D^*^°7r+, D^*^°p+, D^*^\t and ;D^*^°D(*)+ where p+ is reconstructed in tt+ttO mode (|M^+^o - 
Mp+ \ < 0.3 GeV/c^) and a| is reconstructed as p°7r+ (| AfpO^+ - 1 < 0.25 GeV/c^). 

candidates are reconstructed as — > K^n" , K~^TT~7r^, K^-k^ti^-k~ , K^Ti'n^ , 

K'g'K~ Tx'^ Ti'^ and K~K^ . D* mesons are reconstructed by combining the D candidates 
with a pion or a photon. The invariant mass of D candidates is required to be within 
a ±3 MeV/c^ (for D* t:^) and ±10 MeV/c^ (for D* 7) intervals around the nominal D* 
mass. Z)+ candidates are reconstructed in the decay modes ^s^^ ^^"^ K^K~7r~^. 

The invariant mass of the D'^ candidates is required to be within ±15 MeV/c^ interval 
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FIG. 1: Distribution of the beam energy constrained mass in data for fully reconstructed B mesons 
(histogram). The solid curve shows the result of the fit and the dotted one is the background 
component. 

around the nominal mass. D*~^ candidates are defined as -Dj"7 combinations where the 
Z)^7 invariant mass lies in the interval ±15 MeV/c^ around the nominal D*~^ mass. Charged 
B pair events are produced from T(4S') resonance {y/s ~ 10.58 GeV), where the or B~ 
is produced with specific momentum and energy. Selection of the fully reconstructed B 
candidates is made according to the values of two variables : the beam-constrained mass 

Mbc = \/ -^beam ~ Pb the energy difference AE = Eb — Ebeam- Here, Eb and pb are the 
reconstructed energy and momentum of the fully reconstructed B candidate in the center- 
of-mass (CM) system, and -Ebeam is the beam energy in the CM frame. The signal region 
for tagging B candidates is defined as Mbc > 5.27 GeV/c^ and -80 MeV < AE < 60 MeV. 

The Mhr distribution of reconstructed B candidates is fit with the sum of an Argus 
function [lO| and a Crystal Ball function The Argus function models the continuum and 
combinatorial background whereas the Crystal Ball function models the signal component, 
which peaks at the B mass. The purity is defined as S/ {S + B), where S (B) is the number 
of signal (background) events for Mbc > 5.27 GeV/c^, as determined from a fit. Fig. [T]shows 
the Mbc distribution for all Bj-^c candidates in our data set in the signal AE region. In this 
sample, there are cross-feed effects between charged and neutral B. From the Monte Carlo 
simulation, we estimate the fraction of B^ {B~^) events in the reconstruction of B~^ (i?°) 
to be 0.095 (0.090). Then, we obtain Nb+b- = (4.00 ± 0.24) x 10^ and the purity of 0.55, 
where the uncertainty on Nb+ b- is dominated by systematic errors. 

In the events where a i^j-ec is reconstructed, we search for decays into a r plus a neutrino 
and a K plus two neutrinos. Candidate events are required to have one or three signal-side 
charged track(s) with the total charge being opposite to that of the reconstructed B. The 
T lepton is identified in the following decay channels: fi~uu, t~ — > e'uu, t~ Tr'u, 

T~ — > n'Tr^u, and r~ ti'ti^ti'v. 

We require the charged particles to be identified as leptons, pions or kaons. The event is 
required to have zero net charge and -Eecl less than 0.3 GeV where -Eecl is the remaining 
energy calculated by adding the energy of the photons that are not associated with either 
the Bt-cc or the tx^ candidate from r~ — > 7r~7r°z/ decay. For all modes except r~ ir^ir^u 
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TABLE I: Expected background in the signal region for the different selection modes. 



mode we reject events with vr*^ mesons in the recoil against -Brec- 

We place the following requirements on the momentum of the track(s) in the CM, p.,^- > 
0.8 GeV/c for r~ — >• 7r~z/, p^-jrO > 1.2 GeV/c for r~ 7r~7r°z/, p^-^+^- > 1.4 GeV/c for 
r~ — > Tf^TT^TT^u, and px- > 1-2 GeV/c for B^ K^vv. The event is required to have 
the total missing momentum of the event to be greater than 0.2 GeV/c for all modes except 
leptonic decay modes and the direction of missing momentum to be —0.86 < cos6'^igg < 0.95 
in the CM frame. Further requirements are made on the invariant mass of two or three pions 
|M^^ - Mp\ < 0.15 GeV/c^ and |M^^^ - M^+ \ < 0.2 GeV/c^. The selection efficiencies for 
each decay mode we consider in this analysis are determined from a large sample of GEANT- 
based Monte Carlo simulations 3| for B~ t~u and B~ K~ui' events generated by 
EvtGen decay package [l3l| . We compute the efficiency as the ratio of the number of events 
surviving each of our selections over the number of fully reconstructed B"^. 

The most powerful variable for separating signal and background is the remaining energy 
-E'ecl- We use different energy cuts for neutral clusters contributing to the -Eecl for barrel 
part and end-cap parts since the effect of beam background is severe in the end-caps. For 
signal events the neutral clusters contributing to the -Eecl can only come from beam back- 
ground, therefore the signal events peaks at low -Eecl and the background events, which 
contain additional sources of neutral clusters, are distributed toward higher -Eecl values. 

The -Eecl < 0.3 GeV region is defined as the signal region and the 0.45 < -Eecl < 1-5 GeV 
region is defined as the sideband region. The -Eecl shape in the MC distribution is used 
to extrapolate the sideband data to the signal region. The number of MC events in signal 
region and sideband are counted and their ratio {tmc) is obtained. Using the number of 
data in the sideband and the ratio tmci the number of expected background events in the 
signal region is estimated. The background estimation for the different subdecay modes from 
the -Eecl sideband extrapolation is shown in Table |l] The numbers of events in sideband 
region agrees well between MC and data. To obtain the background expected from the MC 
simulation, BB and e^e~ uu, dd, ss, cc events are scaled to equivalent luminosity in 
data. 

The double tag events, for which one of the B mesons is fully reconstructed and the other 
B meson is reconstructed in the set of decay modes B~ D^i~u, where £ is a muon or 
an electron and the is reconstructed in two modes : K^tt~ & K^tt'tt^ti' , are 

used as a control sample to validate the -Eecl simulation. The sources affecting the -Eecl in 
double-tagged events are similar to those affecting the -Eecl distribution in the signal MC 
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FIG. 2: The E-^cL distribution for double-tagged events, plotted for simulation and data. Contin- 
uum background is negligible. 

simulation. The agreement of the -Eecl distribution between data and MC simulation for 
the double-tagged sample, in Fig. |21 is used as a validation of the Eecl simulation in the 
signal MC. 

The main sources of uncertainty we consider in the determination of the B{B^ t^V) 
and B{B~ K~u~V) are the number of B^B~ events with one reconstructed 5, the 
determination of the signal efficiency and the determination of the number of expected 
background events. The number of B^B~ events is determined as the area of the Crystal Ball 
function fitted to the Mbc distribution. Using a Gaussian function as an alternative fitting 
function, we obtain a relative change in the number of events and this difference is assigned 
as the systematic uncertainty on the number of B^B~ events. The main contribution to 
the systematic uncertainties in the determination of the efficiencies comes from uncertainty 
on tracking efficiency, Monte Carlo statistics and particle identification. The uncertainty in 
the expected background comes from Monte Carlo statistics and statistics of sideband data 
events. Estimates of systematic uncertainties are summarized in Table ITTl 

After finalizing the signal selection criteria, the signal region (-Eecl < 0.3 GeV) in the 
on- resonance data is examined. Table 11111 lists the number of observed events in data in the 
signal region, together with the expected number of signal and background events in the 
signal region. Fig. El shows the -Eecl distributions in the data after all selection requirements 
except the one on -Eecl have been applied compared with the expected background. Each 
distribution refers to a different mode. 

Since we do not observe significant excess over the expected background, we set upper 
limits. To extract the upper limits on the branching fraction for B{B~ t~V) and B{B~ — >■ 
K~i>V), we fit the observed Eecl distributions to the expected background and signal, using 
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r decay BR 
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Lepton identification 
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r) r) O O 

A. A I. A 

2.6 
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MC statistics 
Data in sideband 
Total Background Uncertainty 


28.2 31.8 47.5 52.0 57.1 55.8 
12.0 12.2 14.6 27.2 25.0 24.3 
30.6 34.1 49.7 58.7 62.3 60.9 



TABLE II: Systematic uncertainties for the number of B^B^ events with one reconstructed B, 
the determination of the efficiency and the determination of the number of expected background 
events for the different decay channels. 



Decay Mode 
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Efficiency(%) 
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Background Observed 
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3.9 ±0.1 


11.8 ±3.6 
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3.4 ±0.1 
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3.5 ±0.1 
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3.0 ± 1.8 
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1.0 ±0.1 


3.6 ±2.2 
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Total 


33.7 ± 1.4 


13.5 ±0.2 


31.4 ±5.9 


39 




B- 


K vv 




42.8 ± 1.8 


0.70 ±0.03 


2.6 ± 1.6 
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TABLE III: Number of observed data events in the signal region, together with number of expected 
signal and background events. Errors in the background expectation is both statistical and system- 
atic errors. The numbers of expected signal are obtained by assuming that B{B~ — > t~u) = 10~^ 
and B{B- R-vV) = 4 x lO"*'. 
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FIG. 3: ^'ecl distributions in the data after all selection requirements except the one on E-^cL 
have been applied. The red histograms represent the distributions for Monte Carlo simulated signal 
events scaled to B{B- t-V) = 10"^ and B{B- R-vv) = 4 x lO'^. 



maximum likelihood method. The likelihood function C is defined as 



where rib and Us represent the number of background and signal events, respectively, N is 
the number of observed events, iVb is the calculated number of background events, and (Jb is 
the calculated background uncertainty. The variable fs is the normalized signal -Eecl distri- 
bution and fb is the normalized background -Eecl distribution. The negative log likelihood 
function is minimized using MINUIT 0] with respect to Ub for each Ug (= Si ■ Nb+b- • B). 
The 90% C. L. upper limit on the branching fraction B is calculated by 

For — i> T~V, we calculate the likelihood function for each different decay mode {Ci{B)). 
Total likelihood function is defined by 

C{B) = l[Cm (3) 

i=l 

where rich is the number of decay modes for t^V. The full systematic uncertainty must 

be incorporated into the likelihood function. We convolve the systematic uncertainty into 
the likelihood function, C{B), by replacing each point of C{B) by a Gaussian distribution 
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FIG. 4: Likelihood functions for the fit to data after smearing. 

centered at that point with width AB which is determined from systematic uncertainty 
study. To get the value of a particular point of the smeared likelihood function, we integrate 
all the contributions from the Gaussian replaced points of the unsmeared likelihood function. 
To combine likelihood functions of 5 decay modes for t^T^, we simply multiply the 

likelihood functions to produce the combined likelihood. The only complication arises in that 
there are common sources of systematic uncertainty and therefore correlated uncertainties 
between the samples. We throw different normal Gaussian random number for correlated 
and uncorrelated systematics for each systematics sources. In this way, we can smear the 
correlated systematics in the same direction. We obtain the final smeared likelihood function 
by multiplying smeared likelihood functions for B~ t^V. For B~ — > K'vV decay, we 
smear the likelihood function by the total systematic uncertainty. Fig. ^ shows the likelihood 
functions for the fit to data after smearing for B~ — s> t~V (left) and B^ — > K'vV (right). 
We obtain upper limits on the branching fraction at the 90% (C.L.) of 

E{B~ T^V) < 1.8 X 10"^ (4) 

B{B- K-vV) < 3.6 X 10"^ (5) 

In the extension of the Standard Model, one expects significant modification to the B~ — >■ 
T^V decay branching fraction. In the two-Higgs doublet model, the decay can occur via a 
charged Higgs particle. The B~ r~V branching fraction is given as 

B(B- ^ r~T7) = B{B- ^ r-17)sM x r^, (6) 

where th is defined as 

/ tan2/5 2 V (^s 

rH=\l- —^^B , (7) 



m 



H 



rriH is the charged Higgs mass and tan /5 is the ratio of vacuum expectation values of two 



fged 

Higgs doublets |l5|. B{B~ t~I^)sm represents SM contribution given by Eq.((TI). Once 
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FIG. 5: The 90% C.L. exclusion boundaries in the [Mj:/+ , tan /?] plane obtained from the observed 
upper limit on B{B~ t~v). 



we get an upper limit on B{B — > r z/), we can give a constraint on tan/3 and mu- Fig. 
E] shows the 90% C.L. exclusion boundaries in the [mH,tan/9] plane obtained with niB = 
5279 MeV/c^ and B{B- t-v)sm = 0.93 x 10"^ from the CKMfitter prediction compared 



with other experimental searches at LEP[16|, at the Tevatron|17|. and at BABAR. 

In conclusion, we have performed a search for the t~V and B^ K^vV decays 

in a fully reconstructed B sample. Upper limits have been set : 

B{B~ T~V) < 1.8 X 10"^ (90% C.L.) (8) 

B{B- K~vV) < 3.6 X 10"^ (90% C.L.) (9) 

which are the most stringent upper limits on these processes to date. 
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